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NOTICES 


Disclaimers 

The  findings  in  this  report  are  not  to  he  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig¬ 
nated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  he  construed  as  official  Government 
indorsement  or  approval  of  commercial  products  or  services 

referenc'd  herein.  % 


Disposition 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not 
return  it  to  the  originator. 
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ABSTRACT 


The  mathematical  formulation  and  associated  finite-difference  approxi¬ 
mations  for  the  numerical  model  of  a  precipitating  roll  cloud  are  summarized. 
With  that  background,  the  computer  program  used  to  perform  the  computations 
in  experiments  simulating  moist  convect  a  is  described  in  detail. 

The  program  description  includes  disc  ssion  of  the  input  requirements, 
the  intern  .I  operations,  and  the  output.  The  newly  developed  recovery  fea¬ 
ture  of  the  program  is  explained.  All  options,  presently  available  in  the 
program,  are  described.  The  detailed  description  of  the  computer  program  is 
contained  in  individual  flow  diagrams  of  the  main  program  and  the  16  sub¬ 
routines  . 

In  the  concluding  remarks,  the  present  status  of  the  computer  program  is 
given.  Finally,  the  anticipated  future  plans  for  the  program  are  outlined. 
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1.0  INTRODUCTION 


Under  two  consecutive  contracts  with  the  U.S.  Army  Electronics  Command, 
work  has  progressed  for  three  years  toward  the  development  of  a  mathematical 
model  which  simulates  shallow,  moist  convection  culminating  in  a  precipita¬ 
ting  cloud.  The  theoretical  formulation  of  the  model  was  presented  at  the 
outset  of  the  model  development  [1],  Included  in  that  report  is  an  extensive 
dicusssion  of  the  literature  relating  to  convection.  Subsequent  reports  [2, 
3,  4]  have  described  the  successive  developmental  stages  leading  to  the  pre¬ 
sent  model  of  a  precipitating  roll  cloud. 

The  simulation  experiments  are  carried  out  by  numerically  integrating 
the  model  equations  with  suitable  boundary  and  initial  conditions.  The  inte¬ 
grations  are  performed  by  an  e’ectronic  digital  comr-'er.  At  present,  the 
computer  being  used  is  the  UNIVAC  1108  system. 

Because  the  model  has  been  developed  to  the  point  Mint  meaningful  simu¬ 
lation  experiments  can  be  performed,  it  is  felt  that  a  summary  of  the  model, 
and  a  detailed  description  of  the  associated  computer  program  is  desirable. 
This  report  is  primarily  a  description  of  the  rather  complicated  computer 
program  which  controls  numerical  computations  involved  in  the  simulation 
experiments . 

Section  2.0  of  this  report  summarizes  the  mathematical  formulation  of 
the  model.  The  summary  is  composed  of  five  subsections:  2.1  deals  with  the 
peri urbat ion-bas lc-state  approach;  7.2  covers  the  dynamics  ani  thermodynamics 
of  the  mathematical  svstem;  2.3  reviews  the  macrophvsics  associated  with 
water  substance  (vaper,  cloud,  and  precipitation);  2.4  summarizes  the  bound¬ 
ary  and  initial  conditions;  and  2.5  discusses  the  diagnostic  equation  used 
to  determine  the  dynamic  pressure. 
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Section  3.0  presents  a  summary  of  the  finite-difference  approximations 
used  to  perform  the  numerical  computations.  The  finite-difference  approxi¬ 
mations  are  conside.ed  separately  as  they  apply  to  the  interior  points,  to 
the  boundary  points,  and  to  the  initial  time  step.  The  initial  data  is  re¬ 
viewed  in  Section  4.0. 

Section  5.0  contains  a  detailed  program  description.  This  description 
is  presented  in  four  broad  categories.  These  are  input  required  by  the  pro¬ 
gram,  internal  operations  controlled  by  the  program,  output,  and  the  recovery 
feature  of  the  program. 

Section  6.0  contains  concluding  remarks  about  the  program  with  respect 
to  the  present  status  and  future  plans  in  the  utilization  of  the  model. 

Appendix  A  contains  a  description  of  the  symbols  and  Appendix  B  contains 
the  detailed  flow  diagrams  of  the  entire  program. 


2.0  SUMMARY  OF  MATHEMATICAL  FORMULATION  OF  THE  MODEL 


The  theoretical  development  of  the  model  was  treated  in  detail  In  an 
earlier  report  (1).  This  report  will  merely  summarize  the  mathematical  fo-mm- 
lation  to  provide  a  frame  of  reference  for  the  program. 

2.1  The  Perturbation-Basic-State  Approach 

The  shallow  convection,  which  ultimately  leads  to  a  precipitating  roll 
cloud,  is  initiated  by  introducing  a  buoyant  bubble  into  a  model  atmosphere 
initially  at  rest  and  horizontally  uniform.  The  initial  state  of  the  model 
is  considered  to  be  composed  of  two  parts: 

a)  a  basic  state,  and 

b)  superimposed  perturbations. 

The  undisturbed  model  atmosphere  is  the  basic  state.  Therefore,  the  physical 
variables  associated  with  the  basic  state  are  functions  of  the  vertical  coor¬ 
dinate  only.  These  are  potential  temperature  90  ,  density  c0  ,  pressure 
p0  ,  saturation  vapor  pressure  ego  ,  and  relative  humidity  r^  . 

The  initial  superimposed  perturbation  consists  of  the  buoyant  bubble. 

All  departures  from  the  basic  state  become  a  part  of  the  perturbation  fields. 
Therefore,  the  physical  variables  characterizing  the  perturbations  are  func¬ 
tions  of  the  spatial  and  temporal  coordinates. 

The  perturbation-basic-state  approach  provides  a  natural  scaling  of  the 
variables  which  must  be  predicted  by  the  computer.  It  is  well  known  that 
scaling  the  variables  in  a  mult i— iterative  problem  is  quite  useful  in  reduc¬ 
ing  the  possibility  that  round-off  error  will  corrupt  the  solution. 
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2.2  The  Equations  Governing  the  Dynamics  and  Thermodynamics  of  the  Pertur¬ 
bation-Basic-State  System 

The  convective  ceU  is  bounded  vertically  by  two  rigid  horizontal  planes 
a  distance  H  (3  km)  apart.  Laterally,  on  two  opposing  sides,  the  model  is 
bounded  by  two  rigid  vertical  planes  a  distance  2L  (6  km)  apart  and  of  infin¬ 
ite  horizontal  extent.  Midway  between  the  vertical  planes  on  the  lower  of  the 
two  horizontal  planes,  the  origin  of  a  Cartesian  coordinate  system  is  estab¬ 
lished  with  the  x-axis  normal  to  the  vertical  planes  and  the  z-axis  normal  to 
the  horizontal  planes,  as  shown  in  Fig.  1.  All  physical  variables  are  assumed 
always  independent  of  y  ,  and  symmetry  about  the  plane  x  *  0  is  imposed. 

The  symmetry  is  imposed  through  an  initially  symmetric  temperature  perturbation. 
Thereafter,  the  boundary  conditions  and  governing  physical  laws  require  the 
dependent  variables  to  be  either  even  or  odd  functions  with  respect  to  x 
(symmetry  condition).  By  virtue  of  the  symmetry  condition,  the  domain  of  solu¬ 
tion  of  the  numerical  solution  is  0  <,  x  <.  L  ,  0  .<,  z  <.  H  . 

Over  that  domain,  the  dynamic  equations  for  shallow  convection  are  the 
continuity  equation 

V*v  -  0  ,  (2-1) 

and  the  momentum  equation 

-v  D  ^  1  *  .  _  _ . 

—  v  +  vVv  +  —  gk  +  —  Vp  -  0  ,  (2-2) 

dt  P0  Pq 

where  v  is  the  wind  with  components  u  and  w  ,  d  is  the  density,  p 

is  the  pressure,  and  g  is  the  acceleration  of  gravity.  Subscript  zero  indi¬ 
cates  the  basic  state,  and  perturbation  quantities  are  without  subscript.  By 
taking  the  curl  of  Eq.  (2-2),  the  vorticity  equation 
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+  v*Vc 
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(2-3) 
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is  obtained,  where  t ;  is  the  relative  vorticlty  (3w/3x)  -  (3u/3z)  .  Intro- 

^  ^ 

duction  of  tne  stream  function  *  ,  defined  by  v  ■  ui  +  wk  ■  k  *  Vi|i  ,  go 
that 


.  .  ii  w  _  ii 

3z  *  3x 


(2-4) 


satisfies  Eq.  (2-1)  without  further  restriction  on  the  motion  field.  More¬ 
over,  the  vorticity  is  given  by 

4  -  V2*  ,  (2-5) 

where  V2(  )  32(  )/3x2  +  32(  )/3z2  . 

The  stream  function  may  be  predicted  by  means  of  Eqs.  (2-3),  (2-4),  and 
(2-5),  provided  the  density  is  known.  It  can  be  shown  (see  [1],  [3]),  that 
for  shallow,  moist  convection. 
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(2-6) 


In  Eq.  (2-6),  q^  and  q^  are  the  specific  water  contents  in  cloud  and  pre- 
cipitat  ion,  respectively.  It  should  be  noted  that  the  water  substance  vari¬ 
ables  are  the  sums  of  their  respective  basic  states  and  perturbations.  Equa¬ 
tion  (2-3)  can  now  be  rewritten  as 
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(2-7) 


The  potential  temperature  perturbation  is  predicted  by  the  following  thermo¬ 
dynamic  energy  equation 
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(2-8) 
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where  is  the  specific  heat  at  constant  pressure,  where  Qj  is  the  latent 

heat  gained  or  lost  when  water  substance  changes  phase,  and  where  Q2  is  the 
loss  of  heat  due  to  evaporation  of  rain  falling  through  unsaturated  air.  These 
sources  will  be  detailed  in  2.3. 

Given  the  temporal  history  of  the  water  substance  fields,  Eqs.  (2-4), 
(2-5),  (2-7),  and  (2-8)  with  proper  boundary  and  initial  conditions  will  per¬ 
mit  prediction  of  the  motion  and  potential  temperature  fields.  2.4  and  2.5 
will  sumrsarize  the  equations  governing  the  water  substance  fields  and  the 
boundary  and  initial  conditions  for  the  model. 


2.3  Equations  Governing  the  Macrophyslcs  of  Water  Substance 

By  considering  the  continuity  equations  for  water  substance,  it  can 
readily  be  shown  [1]  that  the  equations  for  the  mixing  ratio  of  water  vapor, 
qv  >  the  specific  water  contents  of  cloud  and  precipitation,  q^  and  q^ 
—all  sums  of  basic  state  and  perturbation  quantities — are 


3qv 

dt 
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.  i<k 
Vo  dt 


(2-9) 


39cl 

at 
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+  V*  Vq 
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-  0 


-1 


Vo  dt  +  Vi  +  V? 


Vi  +  V2  +  v3  -  a!  (poqpV)j 


(2-10) 


(2-11) 


where  Sj  ,  S2  ,  and  S3  are  source  terms.  S3  expresses  the  conversion  of 
cloud  to  precipitation;  S2  ,  tl  e  depletion  of  cloud  droplets  by  collision 
with  falling  precipitation  (accretion);  and  S3  ,  the  gain  in  water  vapor 


through  the  evaporation  of  precipitation  falling  through  unsaturated  air.  In 


Eqs.  (2-9),  (2-10),  and  (2-11),  6^  (i  ■  1,  2,  3,  4)  are  symbols  indicating 


t  presence 

<5i- 

1)  or  absence 

(6i  " 

0)  of  certain 

terms . 

The  decision,  as 

whether  each  of 

those  terms  is 

>  present,  is  based 

on  the 

following  criteria: 

61  ‘ 

1 

if 

qcl  >  ° 

or 

q  *  q 
^v  s 

and 

w  >  0 

(2-12) 
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qv  <  qs 
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q  ■  q 

and 

w  <.  0  , 

<5  2  " 

1 

if 

q  >0 

P 

and 

qv  4  qs 

(2-13) 

V" 
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if 

q  >0 

P 

and 

qcl  >  0  . 

63  " 

1 

if 

qcl  "  3 

(2-14) 
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if 

qcl  <  3 
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1 

if 

q  »  0 

P 

and 

qcl  5  0 

(2-15) 

\  “ 

0 

if 

q  »  0 

P 

or 

qcl  •  0  * 

The  sources  in  Eqs.  (2-9),  (2-10),  and  (2-11)  are  given  by 

S j  -  K(cflqcl  -  a)  gm  cm-3  sec”1  ,  (2-16) 

•ita 

S2  -  6.98  *  10~10  En^  coqcl^POqp^  cm_3  sec_1  »  (2-17) 

S,  -  -  1.35  »  10“  1  ‘  c  (q  -  q  )  (c  q  J0*6’  gm  cm”3  sec  1  .  (2-18) 

In  Eqs.  (2-9)  through  (2-18),  V  is  the  relative  fall  speed  of  precipitation, 
a  is  the  threshold  value  of  above  which  a  part  of  the  cloud  is  con¬ 

verted  ro  precipitation  by  autoconversion,  X  is  the  efficiency  of  the  auto¬ 
conversion  process  (autoconversion  parameter),  E  is  the  catch  coefficient 


expressing  the  efficiency  of  the  accretion  process,  nQ  is  the  number  of  pre¬ 
cipitation  particles  per  unit  volume  per  unit  diameter,  and  qg  is  the  satura¬ 
tion  mixing  ratio. 

Returning  now  to  the  heat  sources  in  Eq.  (2-8),  the  source  due  to  phase 
change  is 

dqs 

•  “  Ldt  <2“19) 


and  the  sink  due  to  evaporation  of  precipitation  is  expressed  by 

dqp 

Q2  ’  Ld^  * 


(2-20) 


From  Eq.  (2-20),  it  can  be  seen  that  evaporation  of  falling  precipitation  will 
cause  a  negative  dq^/dt  leading  to  a  decrease  in  0  .  The  mode  of  operation 
of  Qj  is  more  clearly  seen  when  dqg/dt  is  expressed  a9  in  [1], 


where 
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0.622L  -  c  T„ 
_ 2-1 


lcpto  +  j;  «s0. 


*s0 


(2-21) 


(2-22) 


R  is  the  gas  constant  for  moist  air,  Ry  is  the  gas  constant  for  water  vapor, 
L  is  the  latent  heat  of  vaporization,  and  qSQ  is  the  basic-state  value  of 
qg  .  Now  by  Eq.  (2-21),  when  saturated  air  is  moving  upwards,  condensation 
takes  place  and  dqg/dt  is  negative  making  a  source. 

Therefore,  making  use  of  the  relatioiship  already  developed  for  dq^/dt 
with  the  proper  physical  conditions  imposed,  the  thermodynamic  energy  Eq .  (2-8) 
may  be  rewritten  as 
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3t  3z 


S0L  f  dqs  . 

c  T  [41  dt  "  Vo  S 3 

p  0  ' 


(2-23) 


A  further  simplification  in  the  system  of  equations  is  possible  If  Eqs. 
(2-9)  and  (2-10)  are  replaced  by 


If  +  ^‘Vq 


-  o 


-1 


63S1  +  \S2  +  62S3 1  • 


(2-24) 
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+  q 


cl 


qci  -  0  if  q  -  qs 

qci  ”  q  ■  qs  if  q  "  qs 


(2-25) 


Note  that,  with  the  prescribed  conditions,  Eq.  (2-25)  is  a  simple  method  of 
predicting  q^  thereby  eliminating  the  differential  Eq.  (2-10)  for  q  ^  . 
Another  advantage  to  using  Eqs.  (2-24)  and  (2-25)  is  that  in  the  case  of  no 
precipitation,  q  is  conserved,  i.e., 


la  - 

dt 


0 


(2-26) 


The  system  of  equations  which  must  be  solved  for  the  shallow  convection 
experiment  are  (2-4),  (2-5),  (2-7),  (2-11),  (2-12)  through  (2-18)  with  qv 
replaced  by  q  ,  and  (2-21)  through  (2-25). 

2.4  Boundary  and  Initial  Conditions 

The  initial  conditions  are  that  the  model  atmosphere  is  at  rest  and  con¬ 
tains  neither  cloud  nor  precipitation.  A  buoyant  bubble  is  introduced  at  the 
lower  part  of  the  symmetry  boundary.  A  quantitative  description  of  the 
Initial  conditions  will  be  presented  in  Section  4.0  of  this  report. 

The  kinematic  boundary  conditions  follow  from  the  geometry  of  the  model. 
The  rigid  wall  and  symmetry  conditions  require  that  motion  normal  to  the  walls 
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and  the  symmetry  plane  (x  “  0)  vanish.  Also  it  shall  be  required  ‘hat  9p  "  0 
at  the  upper  boundary.  These  conditions  are  shown  in  Fig.  1.  Mathematically, 
the  boundary  conditions  are 


at 

at 

at 


z  *  0  : 

z  -  H  : 

x  -  0  ,  L 


* 

<P 

ip  *  r. 


r,  -  0  ; 

;  -  qp  -  0  ; 
0  . 


(2-27) 


r 


With  these  boundary  conditions,  the  system  of  equations  is  complete. 


2.5  The  Diagnostic  Pressure  Equation 

The  incompressibility  feature  of  the  model  [see  Eq.  (2-1)]  readily  per¬ 
mits  the  derivation  of  a  diagnostic  equation  for  the  dynamic  pressure.  The 
driving  forces  of  the  model  are  the  buoyancy  force  and  the  dynamic  pressure. 
Therefore,  it  is  useful  to  obtain  the  pressure  field  in  order  to  gain  further 
insight  into  tne  dynamics  of  shallow  convection  as  can  be  seen  in  the  latest 
report  [4].  The  departure  of  the  dynamic  pressure  from  the  hydrostatic  pres¬ 
sure  has  been  demonstrated  for  both  the  dry  convection  [2]  and  the  moist  con¬ 
vection  [4]  experiments. 

The  diagnostic  pressure  equation  is  obtained  by  taking  the  divergence  of 
the  equation  of  motion  (2-2)  and  applying  the  continuity  equation  (2-1)  yield¬ 
ing 


'  P 


’cl 


’P 1 


V*(00V7V) 


0 


(2-28) 


i 


i 


The  boundary  conditions  appropriate  to  Eq .  (2-28)  follow  from  application  of 
the  kinematic  boundary  conditions  (2-27)  to  the  equation  of  motion  (2-2). 

Namely,  they  are  that  the  pressure  is  hydrostatic  at  the  upper  and  lower 

* 


11 


boundaries,  and  that  the  pressure  gradients  normal  to  the  lateral  boundaries 
vanish.  Mathematically 


at  z  -  0  ,  H  : 

at  x  •  0  ,  L  : 


bz 

3x 


qd  '.Si 


(2-29) 


Given  the  complete  system  of  equations  summarized  in  the  previous  sec¬ 
tions,  and  the  diagnostic  pressure  equation  with  Itsj  associated  boundary  con¬ 
ditions,  the  shallow  convection  fields  of  motion,  potential  temperature,  water 

f 

substance,  and  dynamic  pressure  may  be  determined  for  any  time. 


3.0  SUMMARY  OF  FINITE-DIFFERENCE  APPROXIMATIONS 

Consider  finite  in.  rements  (Ax,  Az)  on  the  coordinate  axes,  and  let  there 
be  a  mesh  of  points  determined  by  ordered  suras  of  these  finite  increments.  It 
follows  that  that  coordinates  of  any  point  in  the  mesh  are  m/.x  ,  nAz  (m  •  1 , 

2,  ...M,  n  *  1,  2,  Analogously,  discrete  values  of  time  are  considered 

as  ordered  sums  of  a  finite  time  increment  HAt  (i-0,  1,  2,  ...).  Then  the 
system  of  equations  summarized  in  the  preceding  sections  may  be  approximated 
by  finite-difference  equations  defined  on  the  mesh.  The  finite-difference  equa¬ 
tions  may  he  numerically  integrated  with  respect  to  time  between  successive 
time  steps.  The  resulting  system  of  numerically  integrated  equations  may  be 
solved  for  the  dependent  variables  at  any  time  step.  Any  function,  f(x,  z,  t)  , 

l 

defined  on  the  mesh,  will  have  the  value  f  at  the  point  (ra.’.x,  n‘z,  . t )  . 

m  f  n 

Any  basic-state  variable  hQ(z)  ,  being  a  function  of  z  alone,  will  have  the 
value  hr:[n  at  the  point  (nAx,  mAz,  £At). 


3.1  Finite-Difference  Equations  Used  for  the  Interior  Points 

All  of  the  following  equations  are  used  at  the  interior  mesh  points  for 
all  time  steps  except  the  initial  one.  Therefore,  the  domain  of  these  equa¬ 
tions  is  2  <.  m  M-l  ,  2  <  n  i  N-l  ,  i  >.  1  . 

3.1.1  The  centered-dif fercnce  approximations 

The  equations  used  to  approximate  Eq .  (2-A)  to  determine  the  wind 
from  the  stream  function  are 


?. 

u 

m 


^ 

m^n+l  mtn-i 
2Az 


l  l 

"  A1 

!  tn+l.n  m-l.n 

y  m  . — — I..X  A 

m,n  2Ax 


(3-la) 


(3-lb) 


13 


i 


Equation  (2-5)  is  used  to  solve  for  the  stream  function  given  the  vorti- 
city  field  ana  is  approximated  by 


m+1  ,n 


-  2* 


+  «' 


m- 1 .  "i 


(Ax)' 


i  i  i 

it  ..  -  2*  +  *  . 

m  ,n+l _ m.n  m.n-1 

(Az)2 


m.n 


(3-2) 


The  vorticity  Eq.  (2-7)  is  approximated  by  a  modified  centered-difference 
approximation.  The  horizontal  derivative  of  [  (0/9Q)  -  q^  -  q^]  is  approxi¬ 
mated  by  a  vertically  weighted  average  of  centered  differences.  The  equation 
is  numerically  integrated  from  (t  -  l)At  to  (i  +  l)At  .  With  these  modifi¬ 
cations,  the  equation  used  to  predict  the  vorticity  at  (4  +  1 > A t  is 


,i+l 

’m.n 


.1-1  _  A£ 

’m .  n  Ax 


<»,)>  -  (ut.)‘  , 

m-r I  tn  rn—  I  *n 


At  ,  , . t 
—  (w’  ) 

A  z  L  m , 


n+] 


-  («Ol  .1  +  fr1 

m,n-lj  4.x 


-  I  l 

9  —  9 

m+l,n+l  m-1  ,n+l 


0, 


n+1 


+  2 


9*  -  e* 

m-t-l.n  m-l,n 


a  -  f» 

+  m+l.n-1  m-l.n-1 

M°n-1 

.t 


(3-3) 


(\i+l,n+l  ~  ^tn-l.n+l*  ’  2(Ct-l,n  “  ^tn-l.n' 


l  i 

^^m+l,n-l  ^m-l,n-l^ 


l  t  i 

wh  ere  3  "  q  i  +  qD  .An  additional  change  made  in  Eq .  (2-7)  to 

m ,  n  Mn  tn  nn  ,n 

arrive  at  (3-3)  is  the  replaceme.it  of  the  advection  of  vorticity  v*V'.  by 

•* 

the  divergence  of  vorticity  flux  7*'v  .  The  equivalence  of  these  two  terms 
is  olear  in  the  case  of  non-d ivergent  flow. 

The  diagnostic  pressure  Eq.  (2-28)  is  modified  by  replacing  the  last  term 
on  the  left  side  by  its  stream  function  equivalent,  thus 
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V‘0q(v7v)  * 


2p  ffiiL.]2  _iii  lit]  . 
HUx3zJ  ax?  3*2J 


0-4) 


With  this  substitution,  Eq.  (2-28)  is  approximated  by 


p  ..  -  2p  +  p  ,  p  -  2p  +  p  . 

m+l.n  rm.n  rm-l,n  m.n+1  m.n  m.n-1 


2Az  1  ''n+l  6 


m.n+1 


i,n+l  “  C0n-1  6 


m.n-1  _  , l 
eon-i  "  <^»n-1 


+  2o, 


- - -  (ij/  -  2iil  +  i/1  1  (i'1 

n  4(Ax)2  (Az)2  m+l,n+l  m,n  ym- 1 ,  n- 1 '  *■  m-l,n+l 

l  l  1  r  l  l  l  0  S 

-2 ii  +  \l  i-fd.  -  d,1  -  a,  1 

m,n  m+l,n-ly  v  m,n+l  m,n-l  m-l,n  m+l,n 


(3-5) 


The  term  involving  the  stream  function  in  Eq.  (3-5)  is  not  obtained  by  normal 
centered  differences.  It  has  been  found  advantageous  (see  [2])  to  take  cen¬ 
tered  differences  along  the  diagonals  of  the  established  mesh  rather  than  along 
the  axes.  This  procedure  is  valid  only  if  Ax  •  Az  ,  in  which  case,  the  mesh 
length  is  /T  •  Ax  rather  than  Ax  and/or  Az  .  The  advantage  of  using  this 
technique  is  that  it  insures  consistent  truncation  of  both  terms  of  the  right- 
hand  side  of  Eq.  (3-4). 

3.1.2  The  Lax-Wendroff  scheme  applied  to  the  equations  for  heat  and 
water  substance 

It  should  be  noted  that  the  Eqs.  (2-23),  (2-11),  (2-21),  and  (2-24), 

used  to  predict  9  ,  q  ,  q  ,  and  q  ,  respectively,  are  all  of  the  form 

P  s 


+  V-(fv) 


where  f  symbolizes  the  dependent  variable,  and  S  is  the  totality  of  sources 


and  sinks.  These  equations  are  approximated  by  the  two-step  Lax-Wendroff 


scheme 


equation 


For  the  sake  of  brevity,  the  scheme  will  be  described  for  the  general 
above.  For  odd  time  steps,  the  approximation  is 
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till 
f  +  f  +  f  +  f 

m-fl.n  m-l.n  m.n+l  m.n-1 


(3-7) 


In  the  case  of  variables  which  are  retained  as  the  sum  of  the  basic  state  and 
the  perturbation,  e.g.,  q  ,  in  order  to  avoid  smoothing  the  basic  state,  the 
following  expression  is  used  in  place  of  Eq.  (3-7), 

0  f  +  f  ,  +  f ‘  .  +  f  ,  -  2f„  -  f„  -  f~  , 

ra+l.n  m-l,n  m.n+1  m.n-1  °n  °n+l  °n-l 

m,n  4  °n 

(3-8) 

For  even  time  steps,  the  approximation  is 
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m ,  n  Ax 
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Az  [_  m,n-l 

“  (fw>m,n+l. 

+  2AfS4 

m,n 

In  application  to  the  particular  equations  in  9  ,  q  ,  q  ,  and  q  ,  f 

p  s 

needs  only  to  be  replaced  by  the  proper  symbol  and  S  by  the  appropriate 
source-sink  terms. 


This  completes  the  finite  difference  approximations  to  the  system  of 
continuous  equations  required  at  the  interior  points  of  the  finite  mesh. 


3.2  Boundary  Conditions 


The  boundary  conditions,  Eq.  (2-27),  for  the  vorticity  and  atreaa  func¬ 
tion  fields  can  be  readily  expressed  in  finite  difference  form.  They  are 
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(3-10) 


Application  of  these  conditions  to  Eq.  (2-4)  leads  to  one  sided  difference 
approximations 
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These  conditions  are  explicitly  applied  to  the  stream  function  Eq.  (3-1)  and 
the  vorticity  Eq.  (3-2). 

The  boundary  conditions,  Eq.  (2-29),  in  finite  difference  form  are 
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(3-12) 


For  heat  and  water  substance,  the  kinematic  boundary  conditions,  Eq.  (2-27). 
must  be  applied  to  the  respective  differential  equations.  The  resulting  differ¬ 
ential  equations  may  be  of  either  of  the  two  following  forms: 
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or 
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(3-14) 


where  Che  subscript  1  Indicates  the  direction  parallel  to  the  boundary  under 
consideration.  Thus,  for  the  lateral  boundaries,  ?  w  ,  £  z  ,  and  for 

the  upper  and  lower  boundaries,  v^  5  u  ,  x^  J  x  . 

Consider  first  the  Lax-Wendroff  scheme  for  Eq.  (3-13).  For  brevity,  we 
will  demonstrate  the  approximations  only  for  the  lower  boundary.  To  determine 
the  approximation  for  any  other  boundary,  one  need  only  substitute  the  proper 
subscripts  and  wind  components.  The  approximation  for  odd  time  steps  is 


_9  At  _o  f  ?  o 

f  .  +  u ,  (f  ,  -  f *  .)  +  AfS  .  ,  (3-15) 
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where 


fl  +  f 4 
tn+Ll  m-1.1 


(3-16) 


and  the  expression  for  u  is  the  same  as  (3-16)  with  f  replaced  by  u  . 

Oil 

For  even  time  steps,  the  approximation  is 


-fc-1  4 t  -1  -l  .  ,  cl 

m,l  Ax  m,l  m-1,1  m+1,1  m,l 


(3-17) 


This  form  is  applied  to  0  . 

Now  consider  the  Lax-Wendroff  scheme  applied  to  Eq.  (3-14).  For  odd  time 
steps,  the  approximation  is 
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*01-1,1  '  (fu)m+l,l  +  ^m,  1  (um+l,l 


(3-18) 


For  even  time  steps,  the  approximation  is 
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This  form  is  applied  to  q  and  qg  ,  and  to  q^  at  the  lateral  boundaries. 

At  the  upper  boundary,  according  to  the  boundary  conditions,  Eq.  (2-29), 
q^  is  zero.  At  the  lower  boundary,  a  special  scheme  is  employed  to  solve  for 
q^  .  The  form  of  the  differential  equation  is  the  same  as  Eq.  (3-13)  where 
-  (3/3z)(Vq  )  is  included  in  S  .  Instead  of  the  two-step  Lax-Wendroff  scheme, 
a  forward  time  difference  scheme  is  used.  One  other  change  is  that  upwind  dif¬ 
ferences  and  averages  are  used  to  approximate  the  advection  term  rather  than 
the  centered  differences  which  are  employed  in  all  other  approximations.  With 

these  changes,  the  approximations  are:  if  u  ,  >  0  ,  then 
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(3-21) 


This  completes  the  incorporation  of  the  boundary  conditions  into  the 
finite  difference  system. 


3.3  Initial  Time  Step  Approximation 

The  finite  difference  form  for  the  vorticity  Eq.  (3-3)  involves  the  pre¬ 
vious  time  step.  This  cannot  be  done  initially  so  the  following  approximation 
is  used  to  compute  the  vorticity  field: 
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The  finite  difference  form,  Eq.  (3-6),  applies  to  the  remaining  prediction 
equations. 

This  concludes  the  summary  of  the  finite  difference  approximations  for 
the  shallow  moist  convection  simulation  model.  The  next  section  describes 
the  initial  data  required  for  the  experiments. 
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4.0  THE  INITIAL  DATA 


The  initial  data  required  for  a  moist  convection  simulation  are  the 
one-dimensional  basic-state  fields  of  potential  temperature  0Q  ,  density, 
pQ  ,  temperature  T0  ,  pressure  pQ  ,  saturation  vapor  pressure  eSQ  ,  and 
relative  humidity  rQ  ,  and  the  two-dimensional  perturbation  field  of  poten¬ 
tial  temperature  corresponding  to  the  buoyant  bubble.  The  fields  of  TQ  , 
pQ  ,  e3Q  ,  and  rQ  are  used  to  compute  the  initial  fields  of  qg  and  qv 
Initially  the  model  is  assumed  to  be  saturated  and  to  have,  very  nearly 
a  moist  adiabatic  stratification.  The  physical  variables  which  characterize 
this  basic  state  are  presented  in  Table  1.  These  are  the  same  values  which 
have  been  given  in  previous  reports  {3,  4]. 

The  buoyant  bubble  is  the  same  as  used  in  all  previous  experiments  [2, 

3,  4].  It  is  shown  in  Fig.  2  as  a  potential  temperature  perturbation. 


Table  1 


Basic-state  distribution  with  height  of  potential  temperature,  tempera¬ 
ture,  density,  pressure,  saturation  vapor  pressure,  and  relative  humidity. 
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297.00 

297.53 
298.06 
298.58 
299.22 
299.80 
300.37 
300.88 
301.48 
301.99 
302.52 
303.05 
303.61 

j  304.16 
|  304.74 

|  305.35 

I  305.96 
)  306.56 

j  307.07 
j  307.72 
1  308.25 

!  308.81 

j  309.49 
'  310.07 

310.68 
311.28 
311.91 

312.54 
313.05 
313.73 
314.26 


T0CK) 

P0(10~3g/cm3) 

!  P0(mb) 

eSo(mb) 

rQ(%) 

297.00 
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100 
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Fig.  2.  Initial  potential  temperature  perturbation  in  VC. 
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5.0  PROCRAM  DESCRIPTION 


This  section  describes  the  current  computer  program  used  for  the  numerical 
simulation  experiments  In  moist,  shallow  convection.  The  most  meaningful  type  of 
program  description  is  in  the  form  of  flow  diagrams  given  in  Appendix  B.  From 
these,  the  details  of  the  logical  flow  within  the  program  may  be  followed. 

5.1  Input 

The  required  input  to  the  program  can  be  divided  into  three  categories: 
control  variables,  program  constants,  and  initial  fields.  The  latter  were  sum¬ 
marized  in  Section  4.0.  The  others  will  be  discussed  below. 

5.1.1  Control  variables 


The  variables,  which  are  used  to  control  certain  internal  operations 
of  the  program,  permit  flexibility  in  the  computations  and  the  resulting  output. 
By  being  able  to  set  these  variables  at  the  outset,  one  is  able  to  significantly 
modify  the  computational  scheme  without  changing  the  existing  program. 

The  control  variables  in  the  present  program  are  given  in  Table  2.  All 
of  the  control  variables  are  punched  on  one  input  card. 

Table  2 


Control  variables  input  to  the  program 


Variable 

Procram  Identification 

Controls 

MX 

Number  of  horizontal  grid  points 

NZ 

Number  of  vertical  grid  points 

NTS 

Number  of  time  steps 

IDTMX 

Maximum  allowable  time  steps 

ITST 

Frequency  of  relaxing  the  dynamic  pressure  field. 
Pressure  relaxation  always  occurs  initially  and  at 

end  of  first  time  step,  i.e.,  if  ITST  =  3,  relaxa- 

tion  occurs  at  time  step  0,  1,  3,  6,  9,  ... 

JTST 

Frequency  of  printing  output.  Operation  equiva¬ 
lent  to  ITST 

IM0 

Dry  (IM0  =  0)  or  moist  (IM0  =  1)  convection  Simu¬ 
la  t  ion 

w 

M 
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5.1.2  Program  constants 


The  program  constants  are  used  in  the  computations.  By  changing 
these  constants,  except  for  CVN1  and  CVN2,  the  physics  or  geometry  of  the 
experiments  are  modified.  The  program  constants  are  given  in  Table  3. 


Table  3 

Program  constants  inpc  to  the  program 


Program  Constant 

Definition 

CL 

Horizontal  extent  of  model 

CH 

Vertical  extent  of  model 

GV 

Viscosity  coefficient.  The  model  has  so  far  been  run 
with  GV  »  0  (non-vlscous) 

CVN1 

Convergence  criterion  for  stream  function  relaxation 

CVN2 

Convergence  criterion  for  dynamic  pressure  relaxation 

RV 

Water  vapor  gas  constant  (Rv) 

EL 

Latent  heat  of  vaporation  (L) 

CP 

Specific  heat  at  constant  pressure  (cp) 

R 

Moist  air  gas  constant  (R) 

ZN0 

Number  of  precipitation  particles  (n^) 

A 

Threshold  value  for  autoconversion  (a) 

ZK 

Cloud  autoconversion  parameter  (K) 

E 

Precipitation  catch  coefficient  (E) 

V 

Relative  precipitation  fall  speed  (V) 

5.1.3  Input  format 

All  of  the  input  is  on  punched  cards  to  be  read  by  the  computer. 
The  format  of  the  input  cards  may  be  found  in  Table  4. 

Card  1  is  a  recovery  control  card.  The  details  of  its  function  and  for¬ 
mat  are  discussed  in  5.4.2.  It  can  be  seen  from  Table  4  that  a  standard  moist 
convection  experiment  reauires  237  cards  of  input  data.  It  should  be  noted 
that  the  program  can  be  used  to  simulate  dry  convection,  i.e.,  convection  in 
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which  all  water  substance  Is  Ignored.  As  shown  In  Table  2,  the  control  vari¬ 
able  for  a  dry  convection  run  Is  set  equal  to  1.  Therefore,  input  card  5 
must  be  altered  accordingly.  Also,  Input  cards  204-217  are  not  part  of  the 

i1 

Input  required  for  a  dry  convection  experiment. 


Table  4 

Input  card  format 


Input  Card 
Nuobct 

-  1  Fortran 

Contents  j  _ 

,  Format 

i  i 

See  5.4.2 

I  - 

2 

Comments  to  identify  the  experiment;  used  for 
output 

i  12A6 

i 

3 

CL,  CH,  OV 

1 

j  6F12.4 

4 

Blank,  CVN1,  CVN2 

;  5F12.4 

5 

N7. ,  MX,  NTS,  IDTMX,  ITST,  JTST,  1M<S 

715 

6-11 

Basic  potential  temperature  field,  ►"'q 

i  6F12  •  A 

f 

12-17 

Basic  density  field, 

Jn 

!  6F12.4 

18-203 

"Bubble"  (potential  temperature  perturbation  field 

5  ) 

m,n 

6F12.4 

204 

RV,  EL,  CP,  R 

4F12.4 

205 

ZN0,  A,  ZK,  E 

4F12.4 

206 

V 

1F12.4 

207-237 

T  -  ,  p.  ,  es  ,  r 
o  •  n  b  n  n 

4F12.4 

5.2  Internal 

Operat ions 

The  internal  operations  of  the  computer  program  have  been  separated  into 
subroutines  each  of  which  performs  a  specialized  set  of  computations  and/or 
operations.  The  main  program  (MAIN)  controls  the  sequence  in  which  the  sub¬ 
routines  are  called  upon,  thereby  directing  the  general  flow  of  the  computa¬ 
tions.  Each  of  these  aspects  of  the  computer  program  are  summarized  below. 


Within  each  subroutine,  a  block  of  operations  or  computations  is 
performed.  Table  5  states  the  purpose  of  each  subroutine. 


Table  b 

Purpose  of  program  subroutines 


Purpose 

PRTIN 

Prints  the  "dry"  input  to  the  moist  or  dry  convection  simu¬ 
lations 

M0ISIN 

Reads  and  prints  "moist"  input;  computes  and  prints  initial 
q  -  and  q-fields 

RELAX  1 

Performs  relaxation  of  stream  function  field;  selectively 
prints  convergence  information 

RELAX2 

Performs  relaxation  of  dynamic  pressure  field  and  prints 
convergence  information 

HYPRES 

Computes  hydrostatic  pressure  (P2) 

CSBR 

Computes  source  function  for  6-prediction 

PRDCTM 

Computes  source  functions  for  and  controls  prediction  of  q  , 
qg,  and  q  P 

PREDIN 

Predicts  values  of  q  ,  q  ,  and  q  for  interior  points 

P  s 

PREDUL 

Predicts  upper  and  lower  boundary  values  of  q  ,  q  ,  and  q 

P  8 

PREDLR 

Predicts  lateral  boundary  values  of  q^,  qg,  and  q 

EXCNG 

1 

Selectively  controls  flow  of  arrays  to  and  from  auxiliary 
storage 

SM00TH  i 

Transfers  fields  from  dummy  array  (DUMA)  to  proper  array 

DUMP 

Prints  source  functions  for  8  and  q 

P 

SUMPK 

1 

Computes  and  prints  space  averaged  potential  and  kinetic 
energy 

SUMTB  | 

Computes  and  prints  space  averaged  potential  temperature 

PRT0UT  j 

Prints  the  fields  of  potential  temperature  perturbation, 
stream  function,  vorticity,  dynamic  and  hydrostatic  pres¬ 
sure,  specific  water  content  of  cloud  and  precipitation, 
saturation  mixing  ratio,  and  total  specific  water  sub¬ 
stance  (q) 

Some  of  these  subroutines  are  called  every  time  step  while  others  are 
selectively  called,  based  on  Input  values  of  certain  control  variables, 
namely  ITST  and  JTST.  The  subroutines  which  are  selectively  called  are 
DUMP,  RELAX2,  HYPRES,  SUMPK,  SUMTB,  and  PRT0UT.  There  are  also  two  sub¬ 
routines  which  are  called  only  at  the  initial  time  step,  namely  PRTIN  and 
M0ISIN. 

5.2.2  General  flow 

This  section  will  merely  summarize  the  general  flow.  The  summary 
will  neglect  the  variations  under  the  several  recovery  options  as  these  will 
be  discussed  in  5.4.  The  details  of  the  general  flow  are  shown  in  the  flow 
diagram  of  the  MAIN  program,  Fig,  A-l  (see  Appendix  B) . 

The  general  flow  of  the  computations  are  summarized  below: 

1.  Read  input  cards  1  through  203. 

2.  Compute  program  constants. 

3.  Call  PRTIN. 

4.  If  moist  convection  -  call  M0ISIN. 

5.  Call  SUMPK,  SUMTB,  RELAX2 ,  HYPRES,  PRT0UT. 

6.  Compute  initial  vorticity  field. 

7.  Call  RELAX 1. 

8.  Compute  wind  fields  (u  and  w)  from  stream  function  field. 

9.  If  moist  convection  -  call  CSBR. 


10.  Compute  potential  temperature  for  odd  time  steps. 

11.  Compute  potential  temperature  for  even  time  steps. 


Performed 

alternately 


12.  If  moist  convection  -  call  PRDCTM. 


13.  If  this  time  step  is  to  be  printed  -  call  SUMPK,  SUMTB,  RELAX2, 
HYPRES,  PRT0UT. 


14.  Compute  and  test  length  of  next  time  step. 
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’5.  Compute  new  vorticity  field. 

16.  If  this  is  not  last  time  step  -  return  to  step  7. 

17.  ST0P. 

As  pointed  out  in  5.1.1,  the  decisions  which  determine  the  specific 
flow  (dry  or  moist  convection,  print  or  non-print  time  step,  last  time  step) 
are  directed  by  the  input  control  variables. 

This  concludes  the  summary  of  the  Internal  operations  of  the  computer 
program.  The  following  sections  will  discuss  the  output. 


5.3  Output 


There  are  three  types  of  output  from  the  program.  In  order  of  Importance, 
they  consist  of: 

1)  the  physical  and  kinematic  fields, 

2)  fields  or  values  which  either  measure  the  efficiency  of  the  compu- 
tions  or  provide  additional  details  of  the  physical  properties  of 
the  model  (printed  by  the  RELAX 1 ,  RELAX 2 ,  DUMP,  SUMPK,  and  SUMTB 
subroutines),  and 

3)  the  input  data  and  the  initially  computed  fields  (printed  by  the 
PRTIN  and  M01SIN  subroutines). 

5.3.1  Form 

The  primary  output  consists  of  successively  printed  fields  of 
9  ,  <p  ,  t  ,  dynamic  and  hydrostatic  pressure,  ,  ‘Ip  »  9S  .  and  <1  • 

Along  the  margins  of  each  field,  the  x  and  z  indices  are  printed. 

The  secondary  output  takes  several  forms.  The  convergence  information 
printed  out  by  RELAX1  and  RELAX2  consists  of  three  columns  of  figures  corre¬ 
sponding  respectively  to  the  number  of  each  iteration,  the  maximum  residual 


of  th, .t  iteration,  and  the  departure  from  the  convergence  criterion.  Thus, 
a  summary  of  each  relaxation  is  printed  until  convergence  or  50  iterations 
have  occurred.  The  output  from  DUMP  consists  of  printouts  of  source  func- 
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tion  fields  for  0  and  .  The  output  from  SUMPK  consists  of  three  num¬ 
bers,  namely  the  space  averaged  potential,  P  ,  and  kinetic,  K  ,  energies 
and  the  difference  P  -  K  .  The  final  secondary  output  comes  from  SUMTB 
and  consists  of  a  single  number,  the  space  averaged  potential  temperature. 

The  initial  output  from  PRTIN  and  M0ISIN  consists  of  an  orderly  summary 
of  the  input  data,  as  well  as  the  initial  fields  of  0  ,  qg  ,  and  q  . 

5.3.2  Controls  and  options  available 

The  controls  and  options  available  for  the  output,  at  present, 
are  centered  in  the  two  control  variables,  ITST  and  JTST.  By  suitable  choice 
of  ITST,  one  can  specify  the  frequency  of  relaxing  the  dynamic  pressure  field 
and  computing  the  hydrostatic  pressure  field  as  well  as  the  printing  wf  both 
pressure  fields.  The  frequency  with  which  all  other  primary  and  secondary 
output  are  printed  is  established  by  the  value  input  for  JTST.  Therefore, 
the  output  options,  presently  available,  consists  of  a  predetermined  fixed 
frequency  of  pressure  field  output  and  a  predetermined  fixed  frjcuency  of 
other  primary  and  secondary  output.  These  two  frequencies  can  be  either  the 
same  or  different.  The  modes  of  operation  of  ITST  and  JTST  were  given  in 
Table  2. 

5.3.3  Anticipated  changes 

There  are  several  changes  in  the  output  which  are  anticipated. 

These  are: 

1)  water  substance  budget, 

2)  changes  in  q^  output  at  the  bottom  of  symmetry  boundary,  and 

3)  more  output  options. 

The  first  change  will  be  made  shortly.  The  water  substance  budget  will 
consist  of  the  space  averages  of  water  vapor  density,  liquid  cloud  density, 


30 


precipitation  density,  and  cumulative  precipitation  reaching  the  lower  bound¬ 
ary.  All  of  these  will  be  computed  at  selected  tirr’  --t  pr  t"  evil.. at-  i-I.c 
degree  of  computational  agreement  with  the  water  substance  continuity  equation. 

Presently,  the  corner  values  of  remain  zero  for  all  time.  Clearly, 

the  value  at  the  bottom  of  the  symmetry  boundary  should  eventually  become  the 

maximum  value  in  the  q  field.  A  scheme  to  determine  q  at  that  point 

P  P 

will  be  developed  in  the  near  future. 

The  third  anticipated  change  is  more  speculative  than  the  others.  Un¬ 
doubtedly,  it  would  be  desirable  if  the  output  options  were  more  flexible 

than  they  are  in  the  present  program.  If  suitable  criteria  can  be  established, 

ll  I 

it  would  be  most  advantageous  if  the  printing  frequency  could  be  internally 
controlled.  The  criteria  would  be  used  to  determine  significant  changes  in 
the  model.  Then,  detailed  output  could  be  obtained  when  significant  changes 
were  occurring,  whereas  a  more  general  output  would  be  printed  at  a  predeter¬ 
mined  frequency.  If  these  criteria  can  be  established,  more  flexible  output 
options  will  be  developed. 

This  completes  the  changes  in  output  which  are  presently  anticipated. 

The  next  section  will  discuss  the  recovery  capability  of  the  program. 


5.4  Recovery  | 

The  recovery  feature  is  the  most  recent  change  that  has  been  made  in 
the  program.  This  feature  permits  resuming  a  completed  convection  experiment 
from  some  pi ^determined  time  step  under  various  values  of  the  physical  parame¬ 
ters.  The  most  obvious  advmtage  in  having  this  option  available  is  the 
capability  of  assessing  the  sensitivity  of  parallel  experiments  to  varying 
physical  conditions. 
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5.4.1  Method 


By  a  suitable  initial  input  card,  a  given  convection  experiment 
can  be  run  so  that  a  recovery  can  be  accomplished  from  specified  time  steps. 
This  is  achieved  by  storing  on  tape  the  arrays  of  all  unknowns,  as  well  as 
the  necessary  program  constants,  for  the  fixed-interval  time  steps  specified 
on  the  initial  input  card. 

To  recover  from  a  previous  convection  simulation,  by  including  a  suit¬ 
able  initial  input  card,  the  recovery  tape  is  searched  until  the  data  for  the 
predetermined  time  step  are  found.  Then  the  following  sequence  occurs. 

1.  The  data  is  read  from  the  recovery  tape  into  the  computer's  internal 
storage. 

2.  The  physical  constants  (input  cards  2  through  5,  204  through  206) 
are  then  read  into  internal  storage. 

3.  The  general  flow  is  started  at  step  7  (see  5.2.2). 

4.  The  experiment  is  run  to  completion. 

The  details  of  this  sequence  may  be  found  in  the  flow  diagram  of  MAIN  (Fig.  A- 
1).  Note  that  step  2  above  is  the  point  at  which  the  physical  conditions 
for  the  recovery  simulation  may  be  changed  from  those  of  the  original  simula¬ 
tion  experiment. 

5.4.2  Available  options  and  associated  input  controls 

All  of  the  recovery  options  are  guided  by  the  values  of  the  con¬ 
trol  variables  on  the  initial  input  card.  There  are  four  recovery  control 
variables,  namely 

INT  -  designation  of  recovery  tape  unit 

NTPS  -  number  of  time  periods  to  be  skipped  for  recovery  tape  input 
IDINC  -  fixed-interval  for  selecting  time  steps  for  recovery  output 

IBTS  -  starting  time  step  for  recovery  input. 


The  FORTRAN  format  of  the  first  input  card  is  613,  9A6;  and  the  contents  are 
INT,  NTPS,  ID1NC,  IBTS,  blank,  blank,  recovery  comments. 

The  four  available  recovery  options  with  the  associated  values  of  the 
control  variables  are  as  follows. 


Option  1  -  No  recovery  input  or  output 

INT  -  IBTS  -  0 

Option  2  -  Store  necessary  data  on  recovery  tape  every  T  time  step 

starting  at  time  step  0  . 

INT  ■  any  number  between  8  and  26 

NTPS  -  0 

TDINC  -  T 

IBTS  -  0 

Option  3  -  Recover  from  previous  experiment  where  IDINC  »  Tc  at  time 

step  TS  and  store  necessary  new  data  on  recovery  tape 
every  T^  time  steps  (Note:  may  be  different  than  TQ) 

INT  «  any  number  between  8  and  26 

NTPS  -  TS/TQ 

IDINC  -  T 

n 

IBTS  -  TS 


Option  4  -  Recover  from  previous  experiment  where  IDINC  ■  T^  at  time 

step  TS  ,  but  do  not  store  any  data  for  future  recovery 

INT  »  any  number  between  8  and  26 

NTPS  =*  TS/T0 

IDINC  >_  NTS  (number  of  time  steps  in  experiment) 

IBTS  -  TS 


The  derails  of  the  mode  of  operation  under  each  of  these  operations  may 
be  found  in  the  MAIN  flow  diagram  (Fig.  A-l) . 
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6.0  CONCLUSION'S 

This  report  summarizes  the  convection  model  which  has  been  developed  at 
The  Travelers  Research  Center,  Inc.  for  the  sponsoring  agency.  The  primary 
purpose  of  the  report  is  to  give  a  detailed  description  of  the  computer  pro¬ 
gram  which  numerically  Integrates  the  modelling  equations.  The  concluding 
remarks  will  discuss  the  present  status  of  the  computer  program  and  the  future 
plans  for  its  use. 

6.1  Present  Status 

The  results  summarized  in  the  last  report  under  this  contract  [4]  were 
obtained  using  the  program  described  in  this  report  on  the  IBM  7094  system. 
Since  that  time,  the  program  has  been  modified  so  that  it  is  fully  operable 
on  the  UNIVAC  1108  system.  From  the  timing  obtained  on  trial  runs,  it  is 
anticipated  that,  to  perform  a  simulation  of  10  minutes  of  moist  convection, 
will  require  approximately  22  minutes  of  computations  on  the  1108  system. 

This  is  approximately  three  times  faster  than  the  IBM  7094  requires  for  the 
same  simulation. 

Of  the  total  internal  storage  available  for  program  instructions,  approxi¬ 
mately  702  remains  available  after  the  present  program  is  stored.  Also, 
approximately  352  of  the  data  storage  remains  available  when  the  present  pro¬ 
gram  is  running.  A  characteristic  feature  of  UNIVAC  1108  is  that  these  two 
types  of  storage  are  interchangeable  only  at  the  sacrifice  of  speed  in  compu¬ 
tations.  If  this  is  done,  a  considerable  increase  in  data  storage  is  possi¬ 
ble,  and  expansion  of  the  present  program  is  under  consideration. 

6.2  Future  Plans 

The  future  plans  for  this  project  may  be  summarized  in  three  separate 

categories : 
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1)  computational  modifications, 

2)  changes  in  basic  state,  and 

3)  experiments. 

Naturally,  there  is  some  interdependence  between  the  categories. 

The  plans  for  computational  modification  are: 

1)  computation  of  a  water  substance  balance, 

2)  computation  of  variable  q  values  at  lower  corner  of  symmetry 

boundary,  15 

3)  correction  of  space  averaged  potential  energy,  and 

4)  computation  of  a  variable  precipitation  fall  speed. 

The  first  three  modifications  have  been  discussed  in  foregoing  sections  of 
this  report.  The  last  modification  represents  a  change  of  the  present  method 
of  using  a  constant  fall  speed.  A  suitable  formula  obtained  elsewhere,  re¬ 
lating  fall  speed  to  appropriate  cloud  physics  parameters,  will  be  used  in 
the  future. 

The  plans  for  changes  in  the  basic  state  are: 

1)  a  non-saturated  lower  layer,  and 

2)  a  very  stable,  non-saturated  upper  layer. 

These  two  changes  are  quite  simple  to  accomplish.  Moreover,  they  produce  a 
basic  state  more  comparable  to  the  climatological  tropical  atmosphere. 

The  anticipated  experiments  include: 

1)  by  using  recovery,  comparing  the  effects  of  suppression  and  non¬ 
suppression  of  precipitation, 

2)  studying  the  effects  of  he  stable  top  layer,  and 

3)  studying  the  effects  of  the  evaporation  of  precipitation  in  the 
non-saturated  layer  below  the  cloud  base. 

The  next  report  will  discuss  the  results  of  the  experiments  completed 

in  the  next  six-month  period. 
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APPENDIX  A 


List  of  Symbols 

With  the  exception  of  n.  ,  subscript  n  on  any  symbol  is  used  in  this 
report  to  designate  tin-  basic  state  of  the  variable  represented  by  the  symbol. 
Superscrint  ?  and  subscripts  m  and  n  on  any  symbol  in  this  report  are 
used  to  denote  the  value  of  the  variable  represented  by  the  sumbol  at  a  point 
in  the  finite  difference  mesh  with  >'  the  time  index,  m  the  x-index,  md 
n  the  z-lndex.  Table  3  gives  the  alternative  symbols  used  in  the  program  for 
certain  of  the  svrbols  given  below. 
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P2 


threshold  value  for  qc^ 
converted  to  precipitation 
to  0.5  ■  10“  gm  cm"  ' 


above  which  a  part  of  cloud  is 
through  autoconversion;  set  equal 
in  the  experiment  reported  in  [4] 


! 

specific  heat  per  unit  mass  at  constant  pressure  (1.004  •  10 
erg  gm"!  deg"1) 


catch  coefficient  for  depletion  of  cloud  by  precipitation 


saturation  vapor  pressure 

acceleration  due  to  gravity  C981  cm  sec"'  ) 

source  function  for  saturation  nixing  ratio  [Eq.  (2-22)]  and 
vertical  extent  of  the  atmospheric  model;  set  equal  to 
3-10  cm 


parameter  in  expression  for  autoconversion  of  cloud  to  rain; 
set  equal  to  10“  -  sec"1  in  the  experiment  reported  in  [4] 

latent  heat  of  evaporation  <j  2 . 500  ■  10 1 erg  gm" 1 ) ,  and  one- 
half  the  horizontal  extent  of  the  atmospheric  model  in  the 
x-direction;  set  equal  to  3  ‘  10  cm 

parameter  for  size  distribution  of  drops  in  prec ipitation;  set 
equal  to  107  n"~* 


dynamic  pressure 


hydrostatic  pressure 
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'  v  cl 


\i  +  ftP 
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specific  water  content  in  cloud 
specific  water  content  of  precipitation 
saturation  mixing  ratio 
mixing  ratio  of  water  vapor 

gas  constant  for  moist  air  (287  *  104  erg  gm-1  deg-1) 

gas  constant  for  water  vapor  (461.5  «  104  erg  gm-1  deg-1) 

relative  humidity 

source  terms  for  water  substance 

temperature 

t  i  me 

horizontal  component  of  the  wind  in  the  x-direction 

vertical  fall  speed  of  precipitation;  set  equal  to  2  *•  10' 
cm  sec-1  in  thi  experiment  reported  in  [4] 

vertical  component  of  the  wind 

horizontal  coordinate 

horizontal  coordinate 

vertical  coordinate 

(i  =  1,  J,  3,  4)  .symbols  indicating  the  presence  ('.  =  1)  or 
absence  (  •  =  '<)  of  certain  terms 

vorticity  for  two-dimensional  flow 

potential  temperature 

dens i ty 

two-dimensional  stream  function 
unit  vector  directed  along  the  positive  x-axis 
unit  vector  directed  upward  along  the  vertical 
two-dimensional  velocity  vector 

two-dimensional  nabla  operator  [1(3/  x)  +  k(.;/’z)] 
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APPENDIX  B 


Flow  Diagrams 

The  flow  diagrams  of  the  MAIN  program  and  the  16  subroutines  of  the  com¬ 
puter  program  for  the  convection  simulation  model  are  contained  herein.  An 
effort  has  been  made  to  include  a  sufficient  amount  of  detail  to  permit  the 
reader  to  follow  the  logical  flow  of  the  program. 

The  symbols  used  in  these  diagrams  are,  for  the  most  part,  those  which 
have  become  standard  flow  chart  symbols.  The  explanation  of  the  symbols  Is 
found  in  Table  6. 

It  should  be  noted  that  subroutine  SUMPK  presently  applied  only  to  dry 
convection.  This  subroutine  is  being  extended  to  include  moist  convection. 
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Table  B-l 


Flow  diagram  symbols 


Fit.  B-l  (Continued) 
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Write:  NZ.  MX,  NTS,  CL,  CI1,  CVN1.  CVN2 


Fig.  B-4  (Continued) 
B-ll 


MOTt:.  CVS?  ojMH'ihril  In  billing  program 


Fig.  B-5.  RELAX2  -  Performs  relaxation  of  p-field;  prints  convergence 
Information. 
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Fig.  B-S  (Continued) 


ter 


Fig.  B-8.  PRDCTM  -  Computes  source  function  for,  and  controls  predic¬ 
tion  of  q„  ,  q_  ,  and  q. 


Fig.  B-8  (Continued) 


B-21 


i  ."MJv-  •***  ^B^***'1 


Recall  DUMA  from  TEMP 


NOTE:  1.  QR(QR0)  is  a  dummy  variable  for  q(qQ),  q^q  Q),  or  q^ 

2.  IADT  is  specified  by  calling  program  to  modify  the  routine  for 


Fig.  B  9.  PREDIN  -  Predicts  values  of  q  ,  q  ,  and  q  for  interior  points. 


PREDUL  -  Predicts  upper  and  lower  boundary  values 


(MM,  NN,  IWAY,  and  NTAPE  are  specified  in  calling  program) 


NOTE:  MKEY  is  specified  in  calling  program 


NOTE:  1.  ARG  is  a  dummy  argument  for  CSB  or  HQP  as  specified  in  calling 
program. 


2.  LARG  is  specified  in  calling  program. 


RETURN 


Fig.  B-14.  DUMP  -  Prints  source  functions  for  8  and  q  . 


a. 


,\’X- :  .  t  _v„ 


Fig.  B-15.  SUMPK  -  Computes  and  prints  space  averaged  potential  and 
kinetic  engergies  and  their  difference. 
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Fig.  B-15  (Continued) 


i  r  + 

bp=  isnM^lsuMN- — i — J 


PMK  =  BP  -  BK 

.  > 

t 

Write:  BP,  BK,  FN 
kinetic.,  and  potentit 

IK  (potential, 
il-kinetic  energies) 

Fig.  B-15  (Continued) 


(IGO  -  Specified  in  Celling  Progrem)  Compute  reel  time 
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The  mathematical  formulation  and  associated  finite-difference  approximations  for 
the  numerical  model  of  a  precipitating  roll  cloud  axe  summarized.  With  that  back¬ 
ground,  the  computer  program  used  to  perform  the  computations  in  experiments  simulat¬ 
ing  moist  convection  is  described  in  detail. 

The  program  description  includes  discussion  of  the  input  requirements,  the 
internal  operations,  and  the  output.  The  newly  developed  recovery  feature  of  the 
program  is  explained.  All  options,  presently  available  in  the  program,  are  described. 
The  detailed  description  of  the  computer  program  is  contained  in  individual  flow 
diagrams  of  the  main  program  and  the  16  subroutines. 

In  the  concluding  remarks,  the  present  status  of  the  computer  program  is  given. 
Finally,  the  anticipated  future  plans  for  the  program  are  outlined. 
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